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We analyze a detailed model of a Bose-Einstein condensate trapped in a ring optical resonator 
and contrast its classical and quantum properties to those of a Fabry-Perot geometry. The inclusion 
of two counter-propagating light fields and three matter field modes leads to important differences 
between the two situations. Specifically, we identify an experimentally realizable region where 
the system's behavior differs strongly from that of a BEC in a Fabry-Perot cavity, and also where 
quantum corrections become significant. The classical dynamics are rich, and near bifurcation points 
in the mean-field classical system, the quantum fluctuations have a major impact on the system's 
dynamics. 
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I. INTRODUCTION 

In recent years, there has been an explosion of interest 
in optomechanical systems in which at least one degree of 
freedom is cooled nearly to its quantum ground state [IJ- 
4 . In the top-down approach, the mechanical element 
(often one end-mirror of a Fabry-Perot cavity that is al- 
lowed to oscillate) is initially in thermal equilibrium with 
its surroundings and then is cooled via radiation pres- 
sure. On the other hand, in the bottom-up approach, 
the mechanical portion of the system typically consists 
of ultracold atoms trapped inside a high-Q optical res- 
onator. The ultracold atomic system can be a thermal 
sample [5] , a quantum-degenerate Bose-Einstein conden- 
sate (BEC) @[7], or even a quantum-degenerate gas of 
fermions [S]. In the bottom- up situation the mechani- 
cal oscillator(s) are comprised of collective momentum 
modes of the trapped gas, excited via photon recoil [S]- 
[13] . The dynamics of the collective interaction between 
photons and ultracold atoms have been studied in de- 
tail both theoretically and experimentally in the context 
of superradiance in BECs [TH - fS2] as well as collective 
atomic recoil lasing (CARL) [23l430| . In both these situ- 
ations, the interplay between light fields and atomic mo- 
tion leads to feedback effects. The bottom-up approach 
to optomechanics exploits this interplay as well, with the 
added feature that the motion of the atoms is analogous 
to that of a mechanical element driven by radiation pres- 
sure. 

In the case of a high-Q Fabry-Perot cavity the intra- 
cavity standing-wave field couples the macroscopically 
occupied zero-momentum component of the BEC to a 
symmetric superposition of the states with center-of- 
mass momentum ±2hk via virtual electric dipole tran- 
sitions [7] . As discussed in a previous paper [31] , there 
are situations where a ring cavity can lead to atomic dy- 
namics different from the standing-wave situation. This 
is because in contrast to a standing wave, running waves 
permit one in principle to extract "which way" informa- 
tion about the matter- wave diffraction process. As a first 



step toward discussing that question, the earlier work 
considered the difference between classical standing wave 
and counterpropagating light fields, that is, the difference 
in optomechanical properties of condensates trapped in, 
say, a Fabry-Perot and a ring cavity. One main conse- 
quence of the presence of two counterpropagating run- 
ning waves was that in addition to a symmetric "cosine" 
momentum side mode, it becomes possible to excite an 
out-of-phase "sine" mode as well. In the optomechan- 
ics analogy, this indicates that two coupled "condensate 
mirrors" of equal oscillation frequencies but in general 
different masses are driven by the intracavity field. We 
showed that this can result in complex multistable be- 
haviors, including the appearance of isolated branches of 
solutions for appropriate choice of parameters. 

The present paper builds on these results and includes 
two new features. At a classical (operators replaced by 
c-numbers) level, the evolution of the zero-momentum 
mode of the condensate is now also included. Further- 
more, this work also discusses the role of small quantum 
fluctuations in the system, particularly on the occupancy 
of the sine and cosine side modes. As previously dis- 
cussed, together with the original condensate they form 
an effective V-system, with the upper levels - the sine 
and cosine modes - driven by a two-photon process in- 
volving both counterpropagating light fields. At the clas- 
sical level one or the other of these modes can become a 
dark state, but quantum fluctuations will normally pre- 
vent these modes from becoming perfectly dark [32] . It 
follows that measuring correlation functions of the optical 
field provides a direct means to probe the quantum prop- 
erties of the matter-wave side modes. These and other 
aspects of the role of quantum fluctuations are examined 
in the following sections, which consider the situation 
where these fluctuations are feeble and their effect can 
be treated in the framework of a linearized theory. 

This paper is organized as follows: Section II intro- 
duces our model of a quantum-degenerate atomic sys- 
tem interacting with two quantized counter-propagating 
field modes in a high-Q ring resonator, and casts it in a 
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form that emphasizes the optomechanical nature of the 
problem. Section III derives the resulting Heisenberg- 
Langevin equations of motion for the system. It solves 
them first in steady state for the case of classical fields, 
recovering in a slightly different form some key results of 
Ref. [31]. Section IV turns to the study of the role of 
quantum fluctuations on the system dynamics. It starts 
by deriving the linearized equations of motion that gov- 
ern these fluctuations, and then briefly outlines a general 
treatment of quantum correlations applicable in cases 
where many separate noise sources are present. The re- 
mainder of the section exploits these results to analyze 
the different behaviors produced by the interplay between 
classical mean field dynamics and quantum fluctuations 
for selected system parameters. Finally, Section V is a 
summary and conclusion. 



states of momenta 2phk, where p is an integer. For fee- 
ble intracavity fields and large detunings it is sufficient to 
consider the first two momentum side modes, p — ±1. It 
is then convenient to decompose the atomic Schrodinger 
field in terms of its momentum ground state and two 
nearest momentum side modes in terms of the parity 
(rather than momentum) eigenstates 



ip(x) 



CO 

V2 



+ c c cos(2fcir) + c s sin(2fcir) 



. (4) 



where Co, c c and c s are bosonic annihilation operators 
for the zero-momentum component and for the sine and 
cosine side modes of the quantum-degenerate atomic sys- 
tem, respectively. 

To couch the problem in a more transparently optome- 
chanical form, we further make the substitutions 1371 



II. MODEL 



CO 



N 



X + iP 



(5) 



We consider a Bose-Einstein condensate of N two-state 
atoms with transition frequency tu a and mass m, assumed 
to be at zero temperature, confined along the path of 
two counterpropagating optical beams in a ring cavity of 
natural frequency tu c and wave number k c = lu c /c. These 
fields are driven by external lasers of intensity rji and fre- 
quency Lo p . We assume that the atomic transition is far 
off-resonance from the field frequency, so that the upper 
electronic level can be eliminated adiabatically. Neglect- 
ing two-body collisions and in a frame rotating at the 
pump frequency uj c the Hamiltonian for this system is 
then 



H = H. 



opt 



BEC 



(1) 



where 
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H B ec = J dxft(x) ^™ 



h^d^ 

2m dx 2 



tp(x), 



kilo / dxifi (x){d\d\ -\- d\a,2 



Here A = 



+ ala ie 2tkx +a\d 2 e- 2tkx )^(x). 
lu c is the pump-cavity detuning, 



(2) 



(3) 



is the off-resonant vacuum Rabi frequency and go is the 
usual (resonant) vacuum Rabi frequency. 

The photon recoil associated with the virtual tran- 
sitions between the lower and upper atomic electronic 
states results in the population of atomic center-of-mass 
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X\ t 2 + iP\ : 2 



(6) 



The first of these equations is indicative of the fact that 
we assume that the zero-momentum component of the 
atomic sample comprises a macroscopically populated 
component that we treat in mean-field theory via a classi- 
cal amplitude V^, to which are superimposed quantum 
fluctuations resulting from the coupling to the sine and 
cosine side modes. 

The approximate expansion Q and the definitions ([5|- 
^ result in the alternate form of the Hamiltonian 



H ' - H 'opt 



(7) 



where 
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h-^{L e M c + L M s ), 
V 2 



X1X2 + P1P2, 
X1P2 — X2P1, 

VmJd + x Xi + p Pi 



s li 



(8) 



(9) 



where / = {c, s}. In Eq. ^ a constant term has been 
ignored and the energy shift of the atom-light interaction 
has been absorbed into the optical part of the Hamilto- 
nian (hence the primed terms). 
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The operators L and M are quadratic light and matter 
operators, respectively, while the e and o subscripts in the 
light operators indicate their parity under interchange of 
the left- and right-moving light fields. Finally, 



(10) 



is an effective detuning that accounts for the mean-field 
Stark shift of the condensate and 



2hk 2 /m 



(11) 



is the recoil frequency associated with the virtual transi- 
tion. 

As already discussed in Ref. [31] the presence of two 
counter-propagating fields in a ring resonator results in 
a situation that is significantly more complex than is the 
case for a high-Q Fabry-Perot cavity. In particular, the 
optomechanical properties of the condensate are now for- 
mally analogous to those of a system of two coupled mov- 
ing mirrors. This can be seen by exploiting the fact that 
due to the large value of N, we can for now neglect the 
nonlinear terms in M r and M,, so that 




-^L o (V2N + X )+£ ps ,(U) 

where i = {1, 2} and j ' = 3 — i. 

The noise sources are assumed uncorrelated for the dif- 
ferent modes of both the matter and light fields. Because 
the damping originates in the a and c operators, it ap- 
pears in both the X and P equations of motion. In the 
case of the light fields, the noise and damping originate 
from cavity loss, vacuum noise, and laser fluctuations, 
while for the matter fields the primary source of noise 
and damping is 3-body collisions with additional nearby 
non-condensed atoms. In addition, the customary factors 
of \f2n and ^7 multiplying the £'s (see e.g. [33]) have 
been absorbed into their definitions. This simplifies later 
results somewhat. 



and hence 



rr/ 



n VN{L e X c + L X S 



(12) 



(13) 



Thus, rather than having a light-matter interaction pro- 
portional to the light field intensity times the position 
of an effective mirror we now have an interaction with 
two "mirrors" |35] of equal mass and effective oscillation 
frequency, but each of which is driven differently due to 
interference effects between the two counterpropagating 
light fields. 



III. MEAN-FIELD STEADY STATE AND 
STABILITY 

The Heisenberg-Langevin equations of motion of the 
system are easily derived from the Hamiltonian Q , com- 
plemented by appropriate quantum noise and damping 
terms. One finds readily 

Xi = -nX t - APi + Rc(r])V2 

+ ^(McPj + (-l)»M,Xj) + i xi , 



Pi = 
X = 

k = 



V2 

-nPi + AXi + Im(?7)V2 

n 



^/2 

-7*0 



(M c X 1 + (-iyM s P J )+^, 

n 



-7-Po 



V2 
V2 



{L e P c + L P S ) + £ x q, 
{L e X c + L a X s ) + S; p o, 



A. Comparison to previous results 

Before undertaking an analysis of the role of quan- 
tum fluctuations, we investigate the effects of including 
the zero-momentum mode itself as a dynamical compo- 
nent in the classical steady-state of the system. Surpris- 
ingly, perhaps, we find that its inclusion can result in 
the appearance of new dynamical features, such as e.g. a 
Hopf bifurcation. To proceed we neglect all noise opera- 
tors, treat the remaining fields classically in a standard 
mean-field approach, and simply look at the existence 
and stability of the fixed points of the classical version 
of Eqs. (14) as the various parameters are varied. The 



results of one such calculation are shown in Fig. [T] re- 
vealing the dependence of the intracavity photon number 
\ai\ 2 — ^(Xf + P 2 ) on the detuning A. (We use the word 
"photon" rather loosely in the classical description of this 
section; more accurately it is a dimensionless measure of 
the mean light field intensity.) 

When compared to Fig. 2 of Ref. [31], which uses iden- 
tical parameters, we note that a quite similar bistable be- 
havior is observed, with two degenerate stable branches 
the photon number can reach, as before. Which branch 
is reached is dependent on initial conditions and quan- 
tum fluctuations, and, when one field's intensity is given 
by the lower branch, the other's is given by the upper. 
There is one particularly interesting discrepancy when 
the evolution of the zero-momentum mode's occupancy 
is included, however. For a small range of negative de- 
tunings - around A/n —0.5 for the present example 
- the bistable branches become unstable, due to a pre- 
viously unseen Hopf bifurcation [34] ■ We attribute this 
bifurcation to the inclusion of the dynamics of the zero- 
momentum mode, an aspect that was neglected in the 
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FIG. 1: Mean intracavity photon number of both modes as a 
function of the effective detuning A. Here k = 2n x 1.3 MHz, 
7 = 2tt x 1.3 kHz, f2 = 2tt x 3.1 kHz, o; r = 2tt x 15.2 kHz, 
N — 9000, and 771 = 772 = 0.54k. Stable solutions are solid 
and unstable are dashed. 
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FIG. 2: Mean sine mode "position" and "momentum" as a 
function of A. Other parameters as above. The outer mag- 
nitude branches are the position solutions and the inner the 
momentum. The maximum sine mode occupancy N a reached 
is roughly 20. 



analysis of Ref. [3T]. Bistable periodic cycles are present, 
and as before, which is reached depends primarily on the 
initial state of the system. The amplitudes of oscilla- 
tions observed in these stable cycles are quite small for 
the zero-momentum mode when compared to its mean 
value (which is of the order vN), but for the lightly oc- 
cupied side modes, the amplitudes can be comparable to 
the mean values. This is a case where simply making the 
ansatz cq — » yN suppresses certain dynamical features. 

In order to do a meaningful linearized quantum treat- 
ment, we attempt in this paper to avoid parameter re- 
gions where complicated multistable behavior is evident. 
The results of the next section will assist us in this ef- 
fort. Furthermore, a broader search of the parameter 
space has hinted that there are also experimentally real- 
izable regimes in which the classical dynamics goes be- 
yond multistability and becomes chaotic. We hope to 
return to this topic in later work. 



B. Classical dark state 

We now turn to a concrete example in which quantum 
and classical predictions will be shown to differ qualita- 
tively. Specifically, we consider the case of symmetric 
pumping 771 = 772 = i] of the counterpropagating cavity 
modes. We have shown previously in Ref. |31j and con- 
firm in the following that within a classical (for the light 
field) and mean-field (for the atoms) theory the sine mode 
is a dark state, for all but relatively narrow parameter 
regions in which bistability and spontaneous symmetry 
breaking occurs (Fig. [2]). This can be seen easily by re- 
placing all optical and matter- wave operators by classical 
expectation values, for instance 

^0 -> Xo, (15) 

where the overbar indicates the mean. Similar definitions 
are used for the rest of the linear operators. When a 
quadratic operator is written with an overbar, we refer 
to the products of the classical means of its constituents 
(L e = X X X 2 + P 1 P 2 , N c = (X 2 + P 2 )/2, etc.). 

The steady-state solution in this limit is easily obtained 
by setting all time derivatives equal to zero. Because 
the equations of motion are coupled cubic polynomials, 
there may be multiple real-valued solutions. Also, be- 
cause there are 10 variables, we cannot necessarily expect 
to find analytic solutions. Nevertheless, we find that the 
dark sine mode solution always exists and has a closed 
form. For much, though not all, of parameter space, it 
turns out to be the unique, stable mean field solution. A 
primary point of interest for the remainder of this work 
is to identify conditions such that quantum corrections 
lead to a finite occupancy of this mode, in particular in 
or near parameter regions where other solutions appear. 
Thus, we begin with the ansatz 

X s = P s = 0. (16) 

This allows the equations of motion for the matter field 
operators to be simplified, yielding constraints for the 
light field mean values: 

X l P 2 -X 2 P l = 0, (17) 

where |a£)| 2 is the mean number of intracavity photons 
in either of the light modes, assuming the dark sine mode 
ansatz holds. In other words, it is the projected number 
of photons in the cavity modes assuming that the sine 
mode is empty. For completeness and to make more di- 
rect contact with Ref. [3T], we also solve for the other 
classical steady state values. 

x = -n 2 {~f 2 + n 2 )z, 

X c = —^UJrilZ, 

P c = - 7 r!(7 2 + fl 2 )Z, (18) 



5 



1.0 



1.5 



2.0 



lm(r])/K 
1 



-1 
-2 

-3 



FIG. 3: Imaginary part of the pump intensity rj required to 
produce mean intracavity field an- All parameters are as in 
Fig. [l] except r\ itself and A = —1.0 x k. Note that because 
Re(?7) is linear in cto (Re(r;) = hold), V itself is determined 
uniquely as a function of an ■ 



where 



n = n \a D \ 2 V2 



(19) 



is the off-resonant Rabi frequency for the intracavity pho- 
ton number |arj| 2 and 



Z 



(20) 



We note that if we calculate M c and M s , the quantities 
governing light-cosine mode and light-sine mode inter- 
action strengths, we retrieve X C \2N and 0, which are 
exactly the results obtained when occupancy changes in 
the zero-momentum mode are neglected. This shows that 
allowingfor the evolution of the zero- momentum mode as 
in Eq. (I5j) versus fixing cq — > VN will only yield different 
behaviors for the light and side-mode fields in the spon- 
taneously broken symmetry region of parameter space. 
The remaining four equations are easily solved to give 



A,2 



KRe(r))V2 — (A — n X c VN)lm{r])V2 



k 2 + (A- n Q X c VN) 2 
Klm(r?)v / 2+ (A - n X c VN)Rc{r))V2 



(A - n Q x c VN) 



(21) 



We may therefore eliminate without loss of gener- 
ality (thereby requiring cun = ai ( 2 to be real) by appro- 
priate selection of the real and imaginary parts (equiva- 
lently, the phase and strength) of the pumping 77, namely 



Re (77) = nan, 
Im(r]) = -(A-n X c 



N)a D . 



(22) 



Recall from Eqs. (18)-(20) that X c has a nontrivial im- 
plicit dependence on ajj (proportional to a D for small 
old a- n d proportional to a D 6 for large old). Once we be- 
gin the discussion of the quantum fluctuations, we will 



treat «o as the free parameter rather than 77, because it 
allows the quantum equations of motion to be put forth in 
closed form. That such a substitution is possible with- 
out requiring an inordinately large pump intensity 77 is 
shown in Fig. |3j a plot of the relationship between the 
desired intracavity field a d and the imaginary part of the 
required 77. This replacement can therefore allow old to 
be interpreted as the 'predicted intracavity field strength 
for a given pumping strength under the assumption that 
there are no atoms excited to the sine side mode. 



Under the restrictions (22), a bifurcation diagram of 



the steady state solutions of the mean intracavity field 
strengths \ol^\ as a function of ajj are shown in Fig. [4] 
In general, the system has up to six steady-state solu- 
tions, with up to two being stable. For small old, we 
have OL\ — 0L2 — old, but as old is increased the sys- 
tem undergoes a bifurcation, with ol\ 7^ a.% . This oc- 
curs at old ~ -4 for our parameters. This is followed by 
a small window with no stable solution (between about 
old = -45 and .6), at which point, there is again a stable 
symmetric solution where the on are larger than might 
be expected. Presumably, the symmetry breaking cre- 
ates something of a positive feedback loop: a few atoms 
are excited into the sine mode, which shifts the phase of 
the counterpropagating light fields, causing them to in- 
terfere with each other less effectively, thereby increasing 
the net number of photons in the cavity, causing more 
atoms to be excited. Lastly, for much larger light fields, 
the pumping dominates the atoms' symmetry breaking 
ability and the dark sine mode becomes stable again. 
In addition, the character of the bifurca tion diagrams as 
functions of A will change because of (22 1, that is, be- 
cause we use a complex (rather than purely real as in 
Fig. [T]) 77 whose imaginary part is itself a linear function 




FIG. 4: Mean intracavity field |a,| as a function of q_d. All 
parameters are as in Fig. [3] Stable solutions are solid and 
unstable are dashed. The dark sine mode ansatz indeed pro- 
duces a solution over the entire range. However, it is only 
unique and stable for very feeble light fields. As the inten- 
sity is increased, the solution bifurcates and becomes unsta- 
ble, returning to stability for sufficiently high intracavity field 
strength (though it is accompanied by additional unstable so- 
lutions) . 
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FIG. 5: Mean intracavity photon number \oti\ 2 as a function 
of A. All parameters are as in Fig. [l] except rather than a 
fixed rj, we take an = 0.25. Stable solutions are solid and 
unstable are dashed. 
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FIG. 6: Mean sine mode position and momentum as a func- 
tion of A. All parameters as in Fig. [5] The outer branches 
are position solutions and the inner are momentum. 

of A. These diagrams are shown in Figs.[5]-[7] The multi- 
stable behaviors in particular are richer than before. The 
dark sine mode remains stable over the entire range of A 
(straight solid line \cti\ 2 = .0625 = .25 2 at the bottom 
of Fig. [5]), yet it is augmented by an isola consisting of 
additional steady-state solutions with macroscopic side 
mode occupancies. In particular, N s can reach approxi- 
mately 250, almost 3% of the total number of atoms (see 
Fig. [6]). The stability and symmetry of these solutions 
varies greatly depending on the detuning. We will see 
later that quantum fluctuations may affect tunneling be- 
tween these stable solution branches. We remark that for 
the range of dctunings where the isola has a single stable 
solution (between A rj —2.2k and A sa —2.7k) the light 
mode has an additional symmetric stable solution, but 
the side mode still remains dark. It becomes macroscop- 
ically occupied once the symmetry of the optical isola is 
broken. Lastly, we examine the effect of allowing deple- 
tion of the zero-momentum condensate mode; that is, we 
look at the steady-state values X and Pq. From Eq. ([5]), 
we have c = VN + X /y/2 + iP /V2- We find that 
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FIG. 7: Mean Po as a function of A. All parameters are as 
in Fig. |B| 

X (not plotted here) is negative in sign and quite a bit 
smaller in magnitude than yN ; considering only this fact 
we might wonder if the evolution of the zero-momentum 
mode is really worth examination. However, when the 
system's symmetry is broken and the sine mode acquires 
a finite occupancy, P can actually be fairly significant 
compared to the relevant scale, y/2N rj 134 (it should be 
noted that Xq, though small in magnitude, still reduces 
the real part of Co enough such that |co| 2 < N; the con- 
densate cannot "grow" from this treatment). This result 
strongly suggests that the replacement of c by y/N is not 
always the optimal ansatz to make; here it seems that cq 
should be a dynamical, rather than static, quantity. To 
put this another way, we expect the classical steady-state 
solutions of the equations of motion to be comprised of 
coherent states of the light and matter fields, including 
the zero-momentum mode. When the sine mode attains 
a finite occupancy, its presence can cause a major shift in 
the quantum state of the zero-momentum mode beyond 
just the removal of a few atoms. 



IV. QUANTUM EFFECTS 

A. Linearized quantum equations of motion 

We now turn to a discussion of the impact of small 
quantum fluctuations on the dynamics of the system. 
To proceed, and armed with our detailed results on the 
mean-field behavior in hand, we include quantum fluctu- 
ations in the original equations of motion, assuming that 
these fluctuations remain sufficiently small that their ef- 
fects may adequately be described by linearized equations 
of motion. Specifically we introduce the fluctuations in 
the familiar fashion via small quantum corrections to the 
mean-field solutions, 

X ^X + xo, (23) 
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and so forth, and linearize the Heisenberg-Langevin equa- 
tions of motion in these fluctuations. This process is jus- 
tified as long as the quantum fluctuations are small com- 
pared to the classical means (or in the case of the sine 
mode, with its zero mean, as long as the quadratic terms 
are smaller than the linear ones). But we remark again 
that because of the instabilities demonstrated here and 
in Ref. [3T] this will not always be the case, so some care 
must be taken when selecting values for the parameters. 

This linearization procedure yields the ten coupled op- 
erator equations of motion 



Ap t + il B pj 



-KPi 

Xox c 



- A£i — flsXj 

PcPO ~ PoPc - 



" XcXQ 



X 


= -7*o + 


Qp c - 


1- Pc(ii - 


F* 2 ) + 


£x0) 




PQ 


= -JPo - 




- Xc(xi - 


M 2 ) + 








= -jx c + 


w r p c 


+ flpo + 


Pq{xi 


+ x 2 ) - 




Pc 


= -UPc ~ 


UJ r X c 


— CIxq — 


Xo( x i 


+ x 2 )- 




X s 


= -7*s + 


OJrPs 


+ Po{p2 


-Pi)' 






Ps 


= -IPs - 




- Xo(P2 


-Pi)H 


~ £,PS 7 




where i 


= {1,2}, j 


= 3- 
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(24) 



(V2N + X )X C + P P C 



V2 _ 
n a a D (V2N + X ), 
Xc = £laot D X Cl 
Pa = QocudPo, 

p c = n a D P c . 



(25) 



The first two of Eqs. (24) describe the fluctuations of the 
light field, and the last six the matter-wave fluctuations 
in the zero-momentum component and in the sine and co- 
sine side modes. The terms proportional to fls describe 
Bragg scattering between the two counter-propagating 
optical fields due to the material grating formed by the 
zero-momentum matter wave and the cosine mode. The 
coupling between the light and matter operators is deter- 
mined by the constants XcbXoPo and p c , which act as 
small perturbations that couple the evolution of the four 
light and six matter operators. Note that there is no cou- 
pling dependent on the occupancy of the sine side mode 
at the level of these equations of motion. The coupling 
coefficients involve only the classical, mean optical and 
matter wave fields. Thus, the coupling to the sine side 
mode only occurs indirectly via quantum fluctuations in 
the symmetric driving situation considered here. 

While the exact eigenvalues and eigenvectors of the 
10 x 10 matrix defined by these equations cannot be found 
explicitly in closed form, standard perturbation theory 
states that the eigenvalues associated to the uncoupled 
optical and matter blocks of equations, e.g. for \o,c = 



Po.c — 0, match those of the coupled system up to second 
order in the perturbation. Numerical testing confirms 
that the values indeed remain close. These eigenvalues 
are 



Alls 



-K±i(A + Q B ), 
-K ± i(A - fl B ), 

-7 ± *\Jw + y ( Wr ± v/4ft 2 +^) , 
-7±iu; r . (26) 



The reason for this nomenclature is clear when one con- 
siders the corresponding eigenvectors (normal modes of 
evolution) . While their explicit expressions are not possi- 
ble to write in closed form and are exceedingly unwieldy 
even when only expanded to first order in the pertur- 
bation, a qualitative inspection yields some useful infor- 
mation. The two pairs of eigenvalues \l c and \l s cor- 
respond to light-dominated evolution with a small mix- 
ture of cosine and zero-momentum state matter modes, 
in the first case, and of the sine matter mode, in the 
second case. As one might anticipate, the oscillation fre- 
quencies of these optical modes have shifts of opposite 
sign depending on whether they are coupled to the sym- 
metric or antisymmetric matter grating. The four Xmc 
eigenvalues correspond to normal modes dominated by 
the cosine mode and zero-momentum mode (in unequal 
proportion) and coupled to the light fields. The Am« val- 
ues correspond to sine-dominated normal modes coupled 
to the light fields. The sine-dominated normal modes 
only contain a tiny contribution from the cosine mode 
and zero-momentum modes and vice- versa, a direct con- 
sequence of the dark-state nature of the sine mode at the 
classical, mean-field level. This explains why the sine 
normal mode's oscillation frequency is just cu r : to lowest 
order, it is decoupled from the other normal modes of the 
system. Since 7 is typically the smallest dimensional pa- 
rameter, as one increases the coupling between the light 
and matter fields, for example by increasing ao, the real 
part of Amc can become positive even though the change 
in \Xmc\ is relatively small. This leads to the instability 
in the dark sine mode solution seen in Fig. [4] Because 
k>7 the optical transients die out rapidly and the light 
fields follow adiabatically the matter wave fields, with 
noise- and interaction-dominated fluctuations. 



B. A formal parenthesis 

To proceed further we open a small parenthesis to in- 
troduce a somewhat formal result that will prove use- 
ful in the analysis of higher-order correlation functions 
of the quantized atomic and optical fields, in particular 
the correlations of various orders of the matter-wave and 
light modes, as well as the cross-correlations between the 
matter and light fields. This formal development extends 
much of the machinery familiar from systems with a sin- 



gle damped operator to deal with multiple noise sources 
with distinct characteristics. 

Consider a quantum system described, in the Heisen- 
berg picture, by a set of d operators Ok that comprise a 

rf-dimensional vector operator O. Let O evolve according 
to a linear (or linearized) system of d coupled Heisenberg- 
Langevin equations, 



d 
dt 



o(t) = wo(t)+e(t), 



(27) 



where W is a d x d matrix with c-number coefficients and 

£ are noise operators with mean (i.e., any net input to 
the system has already been absorbed into the equations 
of motion). For convenience we have merged the usual 
factors of \[2n needed to preserve commutation relations 

into £, since they may vary for each operator. 

These equations can be integrated in a straightforward 
fashion, yielding 

6{t) = e wt S(0) + f e w ^i(u)du. (28) 



(29) 
(30) 



For each i, j < d, and at all times s, t > we have 

(0,(0)4- (*)) = o, 



where the first condition is satisfied axiomatically, and 
the second one holds for white noise sources, an approx- 
imation that should be adequate for the system under 
study. The expectation values of the operators alone are 
just 



(O(t))=e wt (O(0)) 



(31) 



since (£(*)) = 0. 

Under these conditions we can obtain the correlation 
matrix 



(d(s)®d(t)) = 



„W T t 



(O(0)®O(0))e- 

min(s,t) 

e W(s-«) Ne W-( t -„) dUj(32) 



which lets us compute all quadratic correlations at all 
times, e.g. (x 1 (s)p c (t)) . 

For the rest of the paper, we shall work in the long-time 
limit, in which the transient behavior of the operators 
has decayed to 0. In physical terms, for the model in 
question, this corresponds to times t ^> I/7, which are 
experimentally accessible for long-lived BECs. In this 
limit, the initial values have decayed to irrelevance, and 
the correlations are simply 

/■min(s,£) 

{6{s)®6{t)) = / e w(s - u) Ne wr(t ~ u W (33) 
Jo 

This is the central result of this section. For a Gaussian 
process we can easily determine higher order correlations 



as well: all three-operator correlators are 0, and the four 
operator correlators are given by 

(6(s) <gi 6(t) 6{u) ® 6{v)) = 

(d(s) ® d{t)) ® (6{u) <g> 6{v)) 

+ (d(s)®d(u))®(d(t)®d(v)) 

+ (d(s)®d(v))®(d(t)®d(u)). (34) 

Now, we apply this technique to our model. The op- 
erators O are the x's and p's, with the coefficients W 



given by (24 1. Lastly, we determine the £ x ,v> an d hence 



N from thefollowing relations: 



(L(s)iUt)) 


= 2kS(s 


-t){Nt h - 


HI), 


(35) 




= 2kS(s 


-t)Nt h , 




(36) 


&i(»)&(*)> 


= 2 7 <5(s 


-t)(Nj h - 


-1), 


(37) 




= 2 1 5{s 


-t)Nj h , 




(38) 



where i — {1,2}, I = {0, c, s}, the N th, s are thermal 
noise occupancies of the baths near the characteristic fre- 
quencies of the system as given by Bose-Einstein statis- 
tics, and all other quadratic noise correlations are 0. The 
noise matrix N for the position and momentum opera- 
tors thus has the form of 5 2 x 2 block matrices, each 
with 2N th + 1 for the on-diagonal entries and ±i for the 
off-diagonal entries with all 4 entries multiplied by 7 or 
k as appropriate. Because we are dealing with optical 
photons and a BEC at a temperature of at most a few 
fiK, going forward we take all N t h — > 0, or, cquivalcntly, 
we take the bath temperatures to be 0. We find that 
typically increasing the Nth's just adds directly to the 
occupancy of the corresponding fields. 



C. Second-order correlations and quantum 
occupancies 

With these formal results at hand, and working with 
a set of parameters that are a combination of those in 
Figs. [4] -m and with A = — 1. x n, — 0.25, we are now 
in a position to explore a few results for the cross operator 
correlations, before looking at the quantum-fluctuation- 
augmented occupancy of the side modes. Further calcu- 
lations not presented here have shown that the results 
obtained for these particular parameter values are fairly 
typical of the monostable regime in which we are inter- 
ested. 

As expected, the fluctuations in the zero-momentum 
and cosine modes are virtually uncorrelated with those 
of the sine mode, ((x,p)o tC (x,p) s ) rs 0, but are slightly 
correlated with each other (e.g. (x x c ) w 0.022 for com- 
parison, XqX c 0.11, so the classical mean-field correla- 
tion is a more significant contribution) . By far the largest 
correlation between distinct matter and/or light fields, 
however, is the one that confirms our intuition, namely, 



9 



N S ,N C 
10 r 



8 
6 
4 

2 





0.0 



0.1 



0.2 



0.3 



0.4 



a D 



an ■ This possibility is corroborated by the observed shift 
to the left in the plot of (N c ) versus N c ; that is, when 
quantum fluctuations are included, the cosine mode be- 
haves as if cud were slightly larger. On the other hand, a 
somewhat different behavior is seen in the latter plot. As 
A decreases from and approaches the bifurcation seen 
above in Figs. [5]-[7j the sine mode initially starts to in- 
crease in occupancy, but then its (and the cosine mode's) 
quantum fluctuations are suppressed as |A| increases fur- 
ther. Nevertheless, as we see below, the variance in the 
sine mode's occupation is so large that the quantum fluc- 
tuations may still influence the character of the system's 
behavior in the case —2.5k ~ A ~ —1.0k. 



FIG. 8: Side mode occupancies (N c ), (N s ) as functions of ao. 
Parameters as in Fig. [3] The cosine mode (solid) has a larger 
occupancy than that of the sine mode (dashed). For reference, 
the classical mean N c is plotted as well (dotted line). 

N S ,N C 



-6 




1 



FIG. 9: Side mode occupancies (N c ), (N s ) as functions of 
A. Parameters as in Fig. [5J but the range of A has been 
extended slightly. The cosine mode (solid) again has a larger 
occupancy than the sine mode (dashed). For reference, the 
classical mean iV c is plotted as well (dotted); X c and P c do 
not depend on A so it is constant. 



the sine mode's fluctuations are very strongly correlated 
to those of the light field (e.g. (xiX s ) s» (— l) l 0.32, i = 
1, 2). This shows that, indeed, the occupation of the sine 
side mode is driven by and interdependent with the fluc- 
tuations in the light fields. 

We also consider the occupancy of the side modes as 
functions of the most easily tunable parameters A and 
an . Keeping in mind that all single operator expectation 
values such as (£,-) decay to 0, we have 



(X 2 ) = X 2 + (x 2 ), 



(39) 



etc. When evaluating these quantities we must be careful 
to avoid those regions in parameter space where the dark 
sine mode steady state solution is unstable, in particular, 
we need an < 0.4. The results are shown in Figs. [8] and 
[9] In the former, we see a noticeable occupation of the 
sine mode before the classical bifurcation. This may be 
sufficient to shift the bifurcation point to a lower value of 



D. Variance in side mode occupancy 

Because the system is coupled to a zero temperature 
bath, we compare the variance 

aj = (N?) - (Ntf (40) 

to that of a bosonic system in thermal equilibrium, in 
which case 



2 

°J,TH 



(Nj) 2 + (Nj), 



(41) 



specifically taking the ratio of (|40j) to (41 1. A value less 
than one indicates sub-thermal statistics, as would be 
the case when there is a significant classical mean and/or 
quantum fluctuations are suppressed. On the other hand, 
a ratio greater than one indicates significant fluctuations 
and a matter or light field driven out of thermal equilib- 
rium. These ratios are computed and plotted in Figs. [10| 
and 
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as functions of ac and A, respectively. In the 
former, for weak mean intracavity light field ajj, both 
modes' statistics are thermal in nature. As the applied 
field is increased, the sine mode is perturbed to slightly 
higher variance, whereas the cosine mode at first exhibits 
less than thermal variance, as the classical mean field 
contribution grows. However, the quantum fluctuations 



0" s,c 2 /c s.cth 2 




FIG. 10: Variance of side mode occupancy compared to ther- 
mal variance as functions of «d . Parameters as Fig. [8] Cosine 
mode is solid and sine mode is dashed. 
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FIG. 11: Variance of side mode occupancy compared to ther- 
mal variance as functions of A. Parameters as Fig. [9] Cosine 
mode is solid and sine mode is dashed. 

then take over and its variance grows quickly as ac ap- 
proaches the bifurcation at a value of roughly 0.4. We 
should however take this result with a grain of salt since it 
is at this point that the quantum fluctuation contribution 
to (N c ) exceeds the mean contribution N c , thus endan- 
gering the validity of the linearized treatment. Still, it 
should of course be expected that increased fluctuations 
in the cosine mode significantly alter the existence and 
stability of steady state solutions in this critical region. 

In the second of these figures we plot the variances as a 
function of the effective detuning A. In this case, except 
for detunings very near 0, the cosine mode is almost com- 
pletely classical, as is the sine mode for sufficiently nega- 
tive values of A. This implies that, if a zero-momentum 
condensate is formed and allowed to evolve for the pa- 
rameters given and a A of less than —2.5k or so, and if it 
reaches the dark sine mode steady state, it is quite likely 
to remain there indefinitely, as quantum fluctuations are 
strongly suppressed. But for less negative values of A, 
the sine mode fluctuations are significant. It may be 
possible that these fluctuations "anticipate" the classical 
bifurcation nearby in parameter space, or even that they 
allow the new stable solutions to appear for larger val- 
ues of the detuning than they would otherwise. To test 



this would likely require simulation of the full nonlinear 
quantum evolution of the system. 



V. DISCUSSION AND CONCLUSIONS 

By analyzing a detailed model including two counter- 
propagating light fields and three matter fields, we are 
able to find a region in parameter space, with experi- 
mentally accessible values, where the system's behavior 
differs significantly from that of a BEC in a Fabry-Perot 
cavity, and also where quantum corrections become sig- 
nificant. The classical dynamics are rich, and near bifur- 
cation points in the mean-field classical system, the quan- 
tum fluctuations also have intriguing properties. They 
appear strong enough to shift or perturb the dynamical 
bifurcation points. 

This system's dynamics are richer than the typical op- 
tomechanical system, and they may be exploited in the 
future to investigate numerous non-classical effects. For 
instance, because of the strong cross-correlation between 
the light's and the sine mode's fluctuations, it should be 
possible, by measuring the output light fields, to moni- 
tor the sine mode as it tunnels between different possi- 
ble steady state branches. This would contrast with the 
tunneling suppression seen in, e.g., Ref. |36| . and it relies 
on the optically driven fluctuations in the matter fields. 
Also, for significantly larger condensates, with N « 10 6 , 
apparently chaotic classical dynamics are observed. It 
may be that this behavior persists even with the mini- 
mal (N th — 0) kicks given by the noise operators £. To 
explore these possibilities, further work will be needed in 
the form of a full nonlinear quantum treatment. 
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